First-principles atomistic Wulff constructions for gold nanoparticles 
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We present a computational study for the equilibrium shape of gold nanoparticles. By linking 
extensive quantum-mechanical calculations, based on Density-Functional Theory (DFT) to Wulff 
construction, we predict equilibrium shapes that are in good agreement with experimental observa- 
tions. We discuss the effect of the interactions between a nanoparticle and the encapsulating material 
on the equilibrium shape. As an example, we calculate adsorption of CO on several different Au(hkl) 
and use the results to explain the experimentally observed shape change of Au nanoparticles. 
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Bulk gold is the noblest of all metals^, as demon- 
strated by delicate gold jewels manufactured several mil- 
lennia ago which are found intact in excavations. On 
the other hand, catalysts that include oxide-supported 
gold nanoparticles were found to efficiently oxidize CO 
at room temperature^^ Au is by far the best such 
catalyslP. Among the key factors that determine the ef- 
ficiency of Au catalysts is the shape of Au nanoparticles, 
in particular the 5- and 6-fold coordinated atoms at its 
corner^HU 

The shape of Au nanoparticles has a key role in every 
aspect of their functionality, from sensing^ and biolabel- 
ing applications^ to plasmoniciP^ and photonics^. In 
optoelectronics, quantum leaps between electronic states 
transform light into electricity and vice versa. The 
probability of such a transition depends on the density 
of states and the dipole matrix elements according to 
Fermi's golden rule. For a given size, both wavefunc- 
tions and energies depend critically on the nanoparticle 
shape. For example, the lowest excitation energy for a 
cubic nanoparticle is 10% higher than that of a spherical 
nanoparticle of the same volume^. 

Gold nanoparticles are often found in their equilibrium 
shape. This is a polyhedron enclosed by faces of various 
(hkl) crystal orientations such that the total surface en- 
ergy, 

/^hkllhkli (1) 
hkl 

is minimum. A^u is the total area of faces parallel to 
(hkl) plane of the crystal and ^hkl is the surface tension, 
i.e. the energy required to create a surface of unit area 
that is parallel to the (hkl) plane of the crystal. In order 
to predict equilibrium shape, one needs calculations of 
surface tensions for many different (hkl). Several such 
calculations exist in the literature either based on em- 
pirical potentials^ or limited to Miller indexes of and 



1^, or using quantum-mechanics for low-index faces and 
empirical models for higher indexes^. An accurate and 
systematic calculation of all high-index Au surfaces is 
missing. 

The equilibrium shape is often found to change upon 
exposure to some interacting environment. As Au 
nanoparticles are used in CO oxidation catalysis, CO gas 
is the ideal candidate to test this idea. Changes to shapes 
of higher sphericity upon exposure to CO gas have been 
observed both experimentally^ and theoretically^. The 
interface tension of a metal in equilibrium with a gas is 
found to depend on the surface tension, the adsorption 
energy and the coverage of adsorbates Eq. Q. In or- 
der to predict the equilibrium shape in an interacting 
environment using the Wulff construction method, it is 
necessary to have a systematic calculation of adsorption 
energies for all relevant (hkl) surfaces. 

The Wulff construction has been used to predict equi- 
librium shapes in a variety of systems. Wulff poly- 
hedra are often employed in observations a nd m od- 
els for nanomaterials including Cu catalyst ^ 19 * 20 !, or 
semiconductors^. In the past decade, Wulff shapes 
employing surface tensions from first-principles calcu- 
lations were used for the successful prediction of the 
shape of nanoparticles, including interactions with their 
environmenP2H2U i n the context of ammonia-synthesis 
catalysis, an ab initio determination of a nanoparticle 
shape was used as a first step in the creation of a vir- 
tual nano-catalysP2HH]. In that work, the Wulff polyhe- 
dron was filled with atoms in order to create a realistic 
nanoparticle. The advantage of this method was that 
it allowed for detailed analysis of the atomic positions, 
making it possible to calculate structural quantities such 
as the number of active sites. This virtual nano-catalyst 
was used in other similar reactions, such as ammonia 
decomposition^ Here, we expand this methodology by 
including all possible (hkl) orientations. Moreover, we 
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take into account changes in shape that may be induced 
by interactions between the nanoparticle and its environ- 
ment. We apply our method to supported gold nanopar- 
ticles, a system of high technological importance. 

The paper is organised as follows: in Section I, we 
briefly review Wulff's theory regarding the equilibrium 
shape. In Section II, we present a systematic calcula- 
tion of the surface tension for every Au(hkl) with Miller 
indexes up to 4. In section III, we use these surface ten- 
sions to create atomistic models for Au nanoparticles of 
sizes up to 70 nm, and analyse their structural proper- 
ties, such as the concentration of active sites. In Section 
IV, we generalize our methodology for nanoparticles in 
interacting environment. We provide a simple formula 
that relates the interface tension to the surface tension 
and adsorption energy. We calculate the minimum ad- 
sorption energy of CO on every Au surface with Miller 
indexes up to 3, and use these results to calculate the 
change in equilibrium shape of Au particles upon expo- 
sure to CO gas. We summarize our results in Section 
V. 



I. THE WULFF CONSTRUCTION 

The concept of "equilibrium shape" was postulated 
by Gibbs in the late 19th century. Under thermody- 
namic equilibrium, a given quantity of matter will at- 
tain a shape that minimizes the total surface energy of 
the system. More than a century ago, mineralogist G. 
Wulff proposed that the shape that minimizes Eq. 
is such that the distance of each face from the center is 
proportional to the surface tension of the respective (hkl) 
surfaced 

dhkl ~ Ihkl- (2) 

One begins the Wulff construction by drawing up a 
plane (for example, (111)) at a distance dm from the 
origin followed by planes parallel to (hkl) at distances 
dhkl = diiijhki/^fiii- The equilibrium shape will be the 
polyhedron enclosed by these planes, having thus the fol- 
lowing properties: 

(a) The shape depends on ratios between surface ten- 
sions, and not their absolute values. 

(b) (hkl) planes with high surface tension (usually 
high-indexed ones) will be drawn at greater dis- 
tances and are therefore less likely to appear in the 
equilibrium shape. 

(c) Being steeper, high-index faces are usually hidden 
behind low-index ones, and tend to occupy smaller 
areas in the equilibrium shape even if jhki is low. 

(d) The extra energy associated with the formation of 
edges between two surfaces is not taken into ac- 
count. 




FIG. 1. (Color online) Typical calculated Au nanoparticles 
for various sizes, d: (a) d = 12.1 nm (b) d — 27.2 nm (c) 
d — > oo. In (a) and (b), step and kink atoms are shown in 
darker color. In (c), different colors correspond to different 
kinds of surfaces, (c) was created using WulffmarP 1 . 



(e) The Wulff polyhedron belongs to the same point 
group as the crystal structure of the material. 

In addition to Wulff construction, there exist other 
methods for the study of nanoparticles. Advances in 
computers allow for the direct simulation of nanoparti- 
cles of large sizes using empirical potentials, as done for 
example by McKenn£p21. In that work, a large number 
of different shapes are tested to find the lowest-energy 
one. The Wulff construction is complimentary to that 
method. Here, we use Wulff construction coupled to first- 
principles calculations of surface tensions. This method 
offers a systematic, easy-to-follow recipe for the construc- 
tion of atomistic models of nanoparticles. 



II. SURFACE TENSION OF GOLD SURFACES 

As we are interested in relatively large Au nanoparti- 
cles, we limit our study to nanoparticles where Au atoms 
far from the surfaces are in the ideal fee lattice. This 
is observed in simulations of large clusters^, although 
small gold clusters may have structures very different 
from fcc^. We begin by calculating the surface tension, 
jhkh °f Au by simulations of periodic (hkl) slabs us- 
ing Density- Functional Theory (DFT). We use the open- 
source Dacapo/ASE suite (https://wiki.fysik.dtu.dk). 
We use a plane wave basis with 340 eV cut-off. The 
core electrons are treated with Vanderbilt non-local 
ultrasoft pseudopotentials^. The Brillouin zone of 
the (lll)-(lxl) surface is modelled by a (10x10x1) 
Monkhorst-Pack grid of fc-points; number of ^-points 
in other surfaces is calculated in proportionality to the 
(111) cell. We use the generalized gradient approxi- 
mation (GGA) Perdew-Wang exchange-correlation func- 
tional PW91 for the clean surfaces and the revised 
Perdew-Burke-Ernzerhof functional RPBE for the sur- 
faces covered with CO since this one gives better adsorp- 
tion energies^. For each set of calculations, we use the 
theoretical lattic constant which is found to be 4.22 A 
for RPBE and 4.18 A for PW91, very close to the ex- 
perimental value of 4.08 A. There is a general trend to 
slightly overestimate the lattice constant of noble metals 
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TABLE I. Ratios of surface tensions of Au in comparison to 
other calculations. The last column contains the adsorption 
energy of CO at low coverage on the same (hkl) surface, and 
the adsorption geometry. b=bridge adsorption, t=on-top ad- 
sorption. . energies calculated from the linear fit shown in 
Fig-i 



using GGA33. We model all (hkl) surfaces of fee Au with 
indexes up to 4. Even for CO-covered nanoparticles, no 
(Akl) surfaces are observed in the Wulff construction; for 
this reason we do not consider (5kl) surfaces in this work. 
Atoms in the top two layers from each side are allowed to 
relax, while subsequent slabs are separated by 12 A of 
vacuum. Slab thickness is chosen independently for each 
(hkl) slab until the surface tension converges within 0.01 
J/m 2 . The surface tension, j^ki is derived from 

E s i ab = NE bu i k + 2A'y hkh (3) 

where N is the number of atoms in the slab, E s i a b is the 
total energy of the slab, E bu ik is the energy per atom in 
bulk Au and A is the area parallel to (hkl). 

The results are summarized in Table |Tj Interestingly, 
the ratio of surface energies of different cells is very close 
to the ratio of the areal density of cleaved bonds^H This 
is another example of the unique chemistry of AvW Au 
atoms have a closed (i-shell and have the least preference 
for directional bonds in the entire periodic table. The 
calculated absolute value for 7m is 0.69 J/m 2 , very close 
to 0.64 J/m 2 reported by Wen and Zhangpl and within 
the same order of magnitude as the values report ed by 
state-of-the-art relativistic all-electron calculation d 15 | 16 l 
The nanoparticle shape depends only on ratios between 
surface energies. As shown in Table |TJ our results for the 
ratios between surface tensions agree with more detailed 
calculations within 5% or less. 



III. AU NANOPARTICLES IN 
NON-INTERACTING ENVIRONMENT 

The Wulff construction for Au is shown in Fig. |ljc) . 
It contains 144 vertices and 86 faces of 5 different kinds: 
(111), (100), (332), (211) and (322) in order of total area. 

To construct atomistic models for nanoparticles, we 
start from a large fee crystal. As (111) has the lowest sur- 
face tension, we begin by choosing the number of (111) 
layers. This determines the distance of (111) plane from 
the center of the nanoparticle, dm, and, consequently, 
the nanoparticle size. For other faces, we use Eq. §, 
with the calculated values of jhkl and cut the crystal at 
the correct distances dhki- We calculate the equation of 
the plane defined by every set of three surface atoms, 
and make sure that only faces consistent with the Wulff 
construction appear on the nanoparticle. We consider 
about 30000 different nanoparticles with diameters rang- 
ing from 1.7 nm to more than 100 nm. 

At small sizes, some faces might not be large enough to 
accommodate a single atom, let alone a unit cell of this. 
Very small nanoparticles expose only (111) and (100) 
faces; in particular, our simulated 459-atom nanoparti- 
cle is identical to the one found from simulations and 
X-ray experiments^. In all cases, the shape resembles 
a truncated octahedron consisting mainly of (111) and 
(100) faces, with their edges decorated by several (hkl) 
faces with indexes up to 3. For diameters up to 16.3 nm, 
we find only (111) and (100) faces; as the nanoparticle 
grows in size, different (hkl) orientations start to appear. 
The thermodynamic limit, shown in [ljc) is reached at 
diameters of the order of 100 nm. 

Geometrical features or typical nanoparticles are 
shown in Table [III The area of nanoparticles is calculated 
analytically using the coordinates of vertices; their vol- 
ume is obtained by numerical integration. By fitting over 
a hundred particles of different diameters, we provide 
scaling relations of various properties with total num- 
ber of atoms, in acco rdance with atom-counting models 
for nanoparticle d 39 * 4 ^ . 



IV. AU NANOPARTICLES IN INTERACTING 
ENVIRONMENT 

The equilibrium shape of nanoparticles that interact 
with their environment can be found by means of a Wulff 
construction based on interfacial tensions, 7hw, between 
Au and its environment instead of surface tensions, "fhki- 
It turns out that the two are related by a simple formula: 

litl = 7fc*+*^, (4) 

Slat 

where 9 is the coverage (number of metal-adsorbate 
bonds over number of surface atoms), A at is total sur- 
face area per metal atom and E a( i s the adsorption energy, 
defined as the excess energy per molecule of the system 
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TABLE II. Characteristic data for typical nanoparticles shown in Fig. 1. Shape (a) is typical for particles up to d = 16.3 nm in 
diameter, shape (b) is typical for larger particles and (c) presents fitted values from over a hundred different particles with d < 
16.3 A. N cornt ; r , N e d ge and N sur f are the total number of atoms at vertices, edges and faces of the nanoparticle, respectively; 
Ntotai is the total number of atoms. A is the total surface area and V the volume of the nanoparticle. (hkl) are the appearing 
surfaces in the shape with the percentage of the total area they occupy. 



compared to isolated Au surface and isolated encapsu- 
lating material. Eq. Q includes implicitly the effects 
of adsorbate-adsorbate interactions, as both the adsorp- 
tion energy and the equilibrium coverage depend on such 
interactions. 

To prove Eq. Q , we use the definitions of r f nt and E a ds 
for a slab of metal in equilibrium with some material X: 

E s iab+x = NE bulk + N ads E x + 2A>yjS!i , (5) 

Eslab+X — E s iab + N a d S Ex + N a d S E a ds- (6) 

In the above equations, E s i a b+x is the total energy of the 
slab+X system, Ex is the total energy per molecule of 
X, and N a ds is the number of bonds between slab and X. 
The latter is related to the coverage, 9 and area per sur- 
face atom, A a u by 9 — N ads /N SU rf and A = A at N sur f. 
Substituting into Eqs. ^ and (|6| and using Eq. ([3| 
yields Eq. Q. 

For a typical system (9 = 0.1, E aa - S =0.5 cV), the sec- 
ond term in Eq.Q is about 0.1 J/m 2 , or 10% of jhkl- 
Change in ratios between various jhki will be of the or- 
der of 1%, resulting in very small change in the equilib- 
rium shape. This explains the similarity of nanoparti- 
cle shapes observed in a wide variety of environments: 
our simulations nicely match experimental observations, 
not only for A u clusters^ but also Au particles on C 
nanotubes^ 21 , on TiO^P an d on CeOj 31 . 

On the other hand, shape can change dramatically for 
very small nanoparticles where bonding on face s mig ht 
be very different from bonding on a large surfac e 44 " 5 ^ or 
when small molecules with high adsorption energy are 
adsorbed. The ideal adsorbate to test this idea is CO. 

We calculate the minimum adsorption energy of CO on 
every Au(hkl) with h,k,l < 3. We consider several differ- 
ent adsorption sites to ensure that the global minimum 
is found; as we are interested in very low CO coverage, 
neighbouring CO molecules maintain a distance of more 
than 4.2 A at all cases. In almost every case, CO binds 
atop the lowest-coordinated Au atom with adsorption en- 
ergy being a linear function of th e coordination number 
of this Au atom, z (Fig. [j^HIIlD We use this linear fit 
to obtain adsorption energies for the nine (Akl) surfaces. 



: Au : C •: O 




Z 



FIG. 2. (Color online) Adsorption energy of CO on Au(hkl) 
as a function of the Au coordination number, z, and its linear 
fit. Some characteristic adsorption geometries are shown. 



Adsorption energies and adsorption sites are shown in 
Table ID 

We use calculated adsorption energies together with 
Eq. Q and obtain the equilibrium shape of Au nanopar- 
ticles at low CO coverage shown in Fig. [3] For rough 
surfaces, 'jhkl will be relatively high, but at the same 
time E a ds will be quite low; this results in a compensa- 
tion effect for the two terms in Eq. Q. As the different 
Ihkl are c l° se to each other, the shape has a much higher 
sphericity^ (98%) than the shape in vac uum (93%), in 
excellent agreement with experiment ;) I I 

Exposure of the nanoparticle to CO gas makes it much 
more reactive. This effect has been ob served in first- 
principles simulations of small Au cluster ^ 18 * 48 ) We find 
that the same happens at larger sizes, although it is more 
prominent for smaller nanoparticles. Assuming that all 
step-edge atoms are active, the active-site density dou- 
bles, increasing from about 200 /imol/g to 400 /imol/g 
for the nanoparticle shown in Fig. [3| 

Gold nanoparticles are usually supported on oxides, 
such as MgO or rutile Ti02- The interaction between the 
nanoparticle and the supporting material will also affect 
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FIG. 3. (Color online) Left: Model of a typical Au nanopar- 
ticle (d = 5.4 nm, ca. 5000 atoms) in weakly interacting 
environment (sphericity=93%, 200 /zmol of active sites per 
g) Right: a same size nanoparticle in equilibrium with low- 
pressure CO gas (sphericity=98%, 400 /imol of active sites 
per g). Step and kink atoms are shown in darker color. 



modynamic equilibrium with their environment. Using 
an atomistic version of the Wulff construction, we gener- 
ate Cartesian positions of nanoparticles, which can then 
be used to analyse structural properties. Our results pro- 
vide insight into large nanoparticles that are of interest to 
catalysis, but are inaccessible by direct atomistic simula- 
tions. The calculated nanoparticles match experimental 
results, including the similarity of shapes in weakly inter- 
acting systems as well as the change towards more spher- 
ical shapes upon exposure to reactive gas. The method is 
easily generalized to other materials, and might be use- 
ful for the improved design of nanomaterials with tailored 
physical and chemical properties. 



its shape. The epitaxial growth will introduce strain in 
the nanoparticle^. More important, the the values of 
IhH for the faces attached to the supporting material will 
be very different. A qualitative picture of this interaction 
has been presented by Lopez et alW^. 
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